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It is analytically shown that the asymptotic correlations in exactly solvable models following a 
quantum quench can behave essentially as thermal correlations provided the entanglement between 
two eigenmodes is sufficiently strong. We provide one example and one counter example of this 
observation. The example illustrates the fact that the thermal correlations arise from initial states 
where the entanglement between the eigenmodes stems from the existence of a large energy gap in the 
initial state. On the other hand, the counter-example shows that when the bi-partite entanglement 
of the eigenmodes stems from interactions that do not open a gap, the correlations at asymptotically 
long times are non-thermal. We also show that the thermal behavior concerns only the asymptotic 
correlation functions, as the difference with an actual thermal ensemble can be observed measuring 
the energy fluctuations of the system. The latter observation implies a breakdown of the fluctuation- 
dissipation theorem. 
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I. INTRODUCTION 



Experiments with ultracold gases [iHSl, and, in par- 
ticular, the ongoing efforts to build quantum emulators 
using ultracold atoms loaded in optical lattices, have 
aroused much interest in understanding the thermaliza- 
tion mechanisms of integrable models [ll la-llSj. The lat- 
ter can be used as simple systems to validate a quantum 
emulator by comparing the outcome of the experiment to 
the exact solution, prior to using the emulator to study 
other, more complex models, for which no exact solutions 
are known. 

However, in order to understand the outcome of a sim- 
ulation of an integrable or exactly solvable model it is 
important to understand the effect of the initial condi- 
tions. Since quantum emulators of ultracold atoms are 
largely isolated systems and its evolution is essentially 
unitary, it becomes necessary to understand the condi- 
tions under which the asymptotic state of the system 
can be described by a standard statistical (i.e. Gibbsian) 
ensemble, or, as it was pointed out recently by Rigol and 
coworkers Q, must be described by a generalized Gibbs 
ensemble (GGE). The latter captures the fact that the ex- 
istence of a non-trivial set of integrals of motion strongly 
constraint the non-equilibrium dynamics of the system. 

In Rcf. [l2|, we showed that the GGE can be analyt- 
ically derived for exactly solvable models and a general 
class of initial states. In particular, we showed that de- 
phasing between different modes makes equal time cor- 
relation functions of both local and non-local operators 



non-ergodic, in the sense that in the thermodynamic 
limit, their infinite time limit only depends on the occu- 
pation of the eigenmodes in the initial state only. Thus, 
the asymptotic values of those correlation functions can 
be equivalcntly obtained by assuming that the correla- 
tions with other eigenmodes produce an effective tem- 
perature. This yields a description of the asymptotic 
correlations that is entirely equivalent to the GGE. 

Nevertheless, it was also noticed in Refs. [3,Q that cer- 
tain kinds of initial states can lead to asymptotic values 
of the observables that are essentially indistinguishable 
from those computed with a standard thermal Gibbs en- 
semble. Other cases of (pre-) thermalization have been 
found in integrable [1J| and non integrable [15l - ll7| sys- 
tems for particular classes of initial conditions. Recently, 
Mitra and Giamarchi [l3| also showed that the adia- 
batic introduction of a non-linearity following a quan- 
tum quench in the Luttinger model, can lead to ther- 
malization as described by the standard Gibbs ensemble. 
The authors of Ref. [l^l emphasized the importance of 
"mode coupling" for thermalization. In this work, we 
show that for certain classes of quenches for which two 
sets of modes are strongly entangled in the initial state, 
the GGE ensemble can be arbitrarily close to the stan- 
dard Gibbs ensemble. Using the methods of Ref. p^ . 
we find a simple instance of the mechanism by which ini- 
tial states can lead to correlations that essentially look 
thermal. As we show below, this happens when the eigen- 
modes of the Hamiltonian that describes the evolution of 
a system following a quantum quench have a certain kind 



of bi-partitc quantum correlations (i.e. entanglement) in 
the initial state. As an application, we find an analytical 
explanation for the numerical observations first reported 
in Ref. [3. This fact could be used as a simple proto- 
col to produce asymptotic correlations that essentially 
look as (rather high temperature) thermal correlations 
in exactly solvable systems. Furthermore, we will show 
below that the effective temperature that characterizes 
the asymptotic thermal correlations can be related to the 
entanglement spectrum of a subset of entangled modes, 
which means that the latter is accessible experimentally 
by measuring the effective temperature that characterizes 
the correlations at long times after the quench. 

Before illustrating the above point with one example 
and one counter example, let us first describe the general 
set up which will be addressed below. Consider a system 
containing two subsystems A and B that are initially cou- 
pled together. For times t < 0, the system is described 
by a Hamiltonian of the form Hq = Ha + Hb + Hab 
where Ha,Hb and Hab are quadratic in some eigen- 
modes {ofc, bk} which carry a quantum number k (which 
forms an continuum in the thermodynamic limit) and can 
be fermionic or bosonic, i.e. 



Ha = y^^eA{k)alak, 

k 

HB^J2^B{k)bibk, 

k 

Hab = ^ ^AB{k) albk + b\ak 



(1) 
(2) 
(3) 



The dispersion relations are assumed such that eA{k) 7^ 
eB{k) for essentially all fc, which is required (see below) 
for dephasing between the two subsystems to occur as 
t — > -foo. We can assume that the system is prepared in 
an initial state in contact with a thermal reservoir at a 
temperature T, i.e. po = Z^^ e~^°'^ (such that Trpo = 
1). For t > 0, the coupling between the two subsystems 
Hab disappears, and the two subsystems evolve unitarity 
and uncoupled, according to the Hamiltonian: 



H = H, 



Hf 



(4) 



The existence of the coupling Hab for all i < im- 
plies that in the initial state, po- there are correlations 
(i.e. bi-partite entanglement) between the eigenmodes, 
i.e. {albk) ^0. 

According to the conjecture of Rigol et al. p, the 
asymptotic state of the system can be described by a 
'generalized' Gibbs ensemble (GGE) density matrix that 
is obtained using the maximum entropy principle tak- 
ing into account that the system dynamics is constrained 
by the existence of the set of integrals of motion given 
by Ia{k) = at(fc)a(fc) and h{k) = b'^{k)b{k). The GGE 
density matrix thus obtained reads: 



PGGE = ^GGE °XP ■ 



E 

fc 



\a{k)a\ak + I3{k)b\bk 



where the Lagrange multipliers are determined by the 
initial conditions, i.e. a{k) = ln[(l ± n°- {k)) / n°- {k)] and 
a{k) = \n[{l±n!'{k))ln^{k% with n°(A:) and n^{k) given 
by ^ and ([TU)) (the + applies to bosonic and the — to 
fermonic modes). 

Alternatively, one can arrive at an equivalent result by 
a different route [i2| . Let us first consider the expansion 
of a local operator in terms of the eigenmodes of H: 






Kfe)tfo, 



(6) 

At asymptotically long times after the quantum quench, 
provided eA{k) ^ ^B{k) and certain conditions of 
smoothness are met, dephasing renders the two-point 
correlation function {O'^ {x,t)0{Q,t)) to the following 
form fil: 



hm {0\x,t)0{Q,t)) = Y, [^ti^)]* ^tm4ak) 

k 

+T.['i'^i=^)T^kmbibk). 



(7) 



Thus, we see that the asymptotic correlations of 0{x) 
depend only on the eigenmode occupation in the initial 
state, a behavior that has been termed non-crgodic in 
Ref. [13 • The above sum over k in Eq. (O allows to define 
a mode-dependent temperature for each mode ^^. In- 
deed, this statement is equivalent to the GGE (cf. EqJSI) 
for a broad class of (Gaussian) initial states (see Ref. [l2| 
and below). Thus, it follows that: 

lim {OHx,t)O{0,t)) = {O^{x)O{0))gge- (8) 

i— J--t-oo 

The above result, valid for local operators, can be com- 
bined with Wick's theorem to show that the asymptotic 
behavior of non-local operators is also described by the 
GGE [i2. 

Alternatively, when the correlations are bi-partite, we 
can regard the effective temperature for the modes in 
the subsystem A as due to their entanglement with 
the modes in the subsystem B (and viceversa). Thus, 
whenever we are dealing with (aj^ak) = Tr poaj^ak or 
{b\,bk) = Trpo^I^fcj we can trace out one of the subsys- 
tems, and write: 



n"(fc) = {a\ak) = Tv pAa\ak = TrpcGE a\ak, (9) 
n\k) ={blbk)^TrpBblbk = TrpGGEbibk, (10) 

where pA ~ Ttb Po and pB — T^t^aPo- Therefore, the 
GGE density matrix can be written as: 



Pgge ^ pa® Pb- 



(11) 



(5) 



We can regard the result in Eq. (|lip as a way to relate 
the density matrix of the GGE ensemble to the reduced 
density matrices of the subsystems A and B. Further- 
more, since both pA and pB are hermitian. it is possible 



to write these objects as follows [l9|, |2C 



Pa{b) 



'A{B) 



(12) 



where we have introduced the entanglement Hamiltonian 
of the A (B) subsystem 'Ha(b)^ which is also a hermitian 
operator. Thus, we see pgge is determined by the total 
entanglement Hamiltonian, Ji = T-La + T~Lb . 

The reduced density matrix describing a subsystem A 
of either a pure state or a thermal mixed state plays an 
important role in quantum information theory applied 
to condensed matter systems [1^. For a pure state, the 
von Neumann entropy Sa = " Tr p^i logj pa measures 
the entanglement between two subsystems A and B. The 
latter can be expressed in terms of the entaglement spec- 
trum of Ha- Recently, the von Neuman entropy and 
entanglement spectrum have become an important tool, 
as they can be used to characterize topological quantum 
phases in various kinds of quantum swtems, such like 
graphene [21 [^topological insulators [22], and quantum 
spin chains 23|. In this context, an important question 
that has been addressed in recent times are the conditions 
under the entanglement Hamiltonian T-La can be propor- 
tional to the subsystem Hamiltonian Ha ■ Some examples 
of this fact been discussed in the literature [22, l24| - |27l | . 
As we show in the example below (see section |lT|) , when 
this happens to be the case in a system like the one de- 
scribed above, we can expect the asymptotic correlations 
after the quantum quench to become essentially thermal. 

Using the methods of Refs. [2^, the entanglement 
Hamiltonian Haib) can be determined for a (Gaus- 



sian) initial state of the form po = Zq 



.-Ho/T 



(PO 



|^o)(^o|/(*o|$o)^liere 
at r = 0). Thus m. 



$o) is the ground state of Hq 



nA = J2^n [(1 ± n''{k))/n-{k)] 4a,, (13) 

k 

nB = Y.^n [(1 ± n'ik))/n''ik)] 6^ 6,., (14) 

k 

which, by comparison with Eq. ([5]), allows us to identify 
the Lagrange multipliers a{k) and /3{k) of the GGE with 
the entanglement spectrum of the subsystems A and B. 

Thus, the entanglement spectrum determines the 
asymptotic state following a quantum quench. Similar 
ideas have been discussed by Qi et. al. for the particular 
case of two-coupled edge states using boundary conformal 
field theory [22] ■ Conversely, provided Lagrange multi- 
pliers a{k) and jSik) could determined experimentally, 
we would be able to access the entanglement spectrum 
and the von Neumann entropy of the subsystems A and 
B. However, in actual experiments it may be difficult to 
obtain the full functional dependence of a{k) and /3(fc). 
Thus, below we shall focus on two cases where the en- 
tanglement spectrum takes a simple form, which may be 
easier to measure experimentally. 

The rest of this work is organized as follows. In the 
following section, using the above results, we provide an 



example of the case in which the asymptotic correlations 
arc essentially thermal, which we show to be a conse- 
quence of the entanglement hamiltonians JiAiB) to be 
proportional to the subsystem Hamiltonian, Haib)- Iii 
section mil we provide a counter-example of the fact that 
bi-partite entanglement does not always lead to ther- 
mal correlations. This counter-example illustrates the 
observation that thermal correlations appear provided 
the bi-partite entanglement arises due from a gap in the 
spectrum of the Hamiltonian that determines the initial 
state. Finally, in section llVi we provide a discussion of 
our results and show that, even when the correlations 
look essentially thermal, there are certain observables 
like the energy fluctuations that still differ from their 
thermal values, a fact that signals a breakdown of the 
fluctuation-dissipation theorem in the asymptotic state. 
The appendix contains some technical details regarding 
a continuum version of the model discussed section [III 



II. EXAMPLE 

Let us first consider a model that has been numerically 
studied earlier by Rigol and coworkers [3] . The model de- 
scribes a system of hard-core bosons in ID that initially 
(i.e. for t < 0) move in the presence of superlattice po- 
tential. The hard-core bosons in ID can be treated using 
a Jordan Wigner transformation [23,l30|, which maps the 
hard-core bosons to non-interacting fermions and, in the 
case of a superlattice of strength A, leads to the following 
quadratic Hamiltonian: 

L L 

where / and /-,- are creation and annihilation operators 
for spinless fermions at site j [j = 1, . . . , L, for a lattice of 
L sites). Rigol et al. |7[ numerically found that, starting 
from the ground state of i^o, if the superlattice term 
oc A is suddenly switched off at t = 0, and the system is 
allowed to evolve unitarity according to 



H 



L 

E 



(/]/, + !+ /j+i/,), (16) 



the long-time behavior of the momentum distribution can 
described by a standard Gibbs canonical ensemble, 



z 



(17) 



for which the effective temperature, Tes, was found to ap- 
proach A/2 for A > 1. In what follows, we will analyti- 
cally demonstrate that this numerical observation indeed 
follows from the existence of a strong bi-partite entangle- 
ment between two sets of eigenmodes of H . 



We begin by Fourier transforming Hq and H by us- 
ingdl 



1 



-ikxi 



Vl 



h. 



(18) 



j=i 



with Xj = j, the Hamihonian (|15|) can be written as 

ifo = ff + A E {flf^+- + fl+J^) ' (19) 

k 

H = Y.^kiflh-fUjk+.) (20) 



where oj^ = — 2cosfc and — 7r/2 < k < tt/2. The Hamil- 
tonian describing the state the system at t < 0, namely 
Hq, can be brought to diagonal form by means of the 
following canonical transformation: 



7fc = cos Ok fk + sin6'fe/fc+^, 
4 = - sin Ok fk +cos9kfk+7v: 
with tan 26k ~ — ■ Hence, 

k 



(21) 
(22) 



(23) 



where Ek = ya^| + A^. Note that the transformation in 
(PTj) and (22) implies the existence of strong bi-partite 
quantum correlations (i.e. entanglement) between the 
modes at k and k + ir, which manifest in, e.g. (flfk+Tr) = 
-isin20fc = -2|:^O. 

As discussed in Refs. [^, ll2|, the asymptotic momen- 
tum distribution of the hardcore bosons for t -^ +oo 
can be obtained from the Fourier transform of the one- 
particle correlation function of the bosons, which in 
turn can be written as a Tocplitz determinant involving 
correlation two-point correlations of the Jordan- Wigner 
fermions: 



lim g^ ' {xi — Xj, t) 

t— >-l-oo 






ao 



a-ri+2 



where 



a-t-j+i 



lim {A,{t)B,{t)). 



i— > + oo 



(24) 



(25) 



The thermodynamic limit is implicitly understood in the 
above expressions; we have also introduced the notations: 



Aj = /J -I- fj and Bj = /J — fj. The above correlation 
function of the (local) operators Ai and Bj can be shown 
to be: 

«^-j+i = ^ E e~''^"'""^-^ [2n(fc) - 1] 

k 



(-1) 



*-J 



L 



^ g-ifc(x.-:r,) J2„,(/, -H tt) - 1] (26) 



Next, we use that (K = Q,tt): 

n{k + K) = {fl+j,fk+K) = Tr PoII+k^+k (27) 

= TrpKfl+Kh+K, (28) 

where po = |*i'o)(*!'o|- 'I'o being the initial state, that is, 
the ground state of Hq (cf. Eg. p^ ) and the reduced 
density matrices: 



Pk — Trfegs_^po, 



(29) 



where 5*0 = (— 7r/2, +7r/2], and S^^ = (— tt, +7r] — Sq. In 
other words, the asymptotic correlations can obtained 
from the reduced density matrix resulting from tracing 
one one the two sets of entangled modes with k belonging 
to either 5*0 or St^- 

As explained in section HI the reduced matrices pk=o,ti 
can be obtained analytically in terms of the occupation 
numbers of n{k + K) = {fl,j^fk+K) in the initial state. 
Using dSI]) and jlSl), we find: 



n{k) 



n{k - 



(30) 
(31) 



2 V Ek^ 

1 / l^k 
- 1 + — 

2 V Ek, 

and, following the discussion in section HI the Lagrange 
multipliers determining the GGE density matrix read: 

a{k) ^In [{I -n{k))/n{k)] 

Ek + LOk 




7r))/n(fc -t- tt)] 



(32) 



For A 3> Wfe, Ek can be approximated by A, and there- 
fore a{k) ~ 2wfe/A ~ -/3(fc). Thus, the GGE density 
matrix, Eq. ([5|) reduces to a standard Gibbs ensemble, 
Eq. dni) with 



Th 



A/2, 



(33) 



which is in agreement with the numerical observations 
of Rigol and coworkers [7|. However, it is important 
to note that the above thermal ensemble and the re- 
sult of Eq. p3p is only an approximation to the actual 
GGE ensemble determined by the Lagrange multipliers 
in Eq. p2[). However, this approximation becomes better 
and better for larger values of A, which implies that nu- 
merically (and experimentally) the GGE and a standard 
thermal Gibbs ensemble will be essentially indistinguish- 
able. 

It is worth noting that the above results are also rele- 
vant for a special limit an integrable field theory in one 
dimension, namely the sine-Gordon model: 



H. 



sG 



Hq '-J- I dx cos 2(j), 



Hq = - I dx 



K-'{d,(l)Y + K{d, 



(34) 
(35) 



where ao is a short-distance cut-off, and the phase and 
density fields 0{x) and (f>{x), are canonically conjugated 
in the sense that the obey: [4>{x),dx'd{x')] = iTrS{x — 
x'); V is the speed of sound and K is a, dimensionless 
parameter that determines the ground correlations of the 
system. Indeed, in equilibrium, the model exhibits two 
phases for g > 0, namely, a gapped phase for i^ < 2 and 
gapless phase ioi K > 2 [231 ■ Foi" K = 1, as described 
in the the Appendix, this model can be mapped onto a 
model of massive fermions in ID: 

p 
-f A ^ [^i{p)Mp) + i'l{p)i'Rip)] , (36) 



with Ep = vp and A oc _g. The model in Eq. p6p can be 
regarded as the continuum limit of the model of Eq. ([T9)) . 
For large values of the gap, A ^ va^^ , the multipliers 
a{p) and /3(p) (see calculation in the Appendix) become 
proportional to the dispersion relation of the Hamiltonian 
that governs the time evolution, ±6° respectively, with 
an effective temperature T^s (it is the same for the two 
branches of fermions) given by 



Teff 



2tanh 



(37) 



2T 



We notice that the effective temperature depends on the 
temperature of the initial thermal state. This effective 
"final" temperature is always higher than the initial tem- 
perature because it contains the bi-partite entanglement 
between k modes. For high initial temperatures the effec- 
tive temperature to which the system thermalizes is the 
same as the initial one. In the case of a zero-temperature 
initial state, the effective temperature results T^s — A/2 
similarly to the case of XY model studied previously. 

Generally speaking, the analysis described above is also 
related to the discussion of the conditions under which 
the entanglement Hamiltonian JiAiB) s-^d the Hamilto- 
nian of the subsystem Hy^rm are (approximately) pro- 
portional to each other. If this is the case, then, ac- 
cording to the discussion of section HI the GGE density 
matrix, Eq. [SJ will be well approximated by the ther- 
mal density matrix of (|17p . Indeed, recently Peschel 
and Chung [27[ addressed the problem of the propor- 
tionality between JiAiB) &iid H^tgy By considering a 
model of two two species of fermions with opposite dis- 
persion Wfc and coupling A, which leads to an energy 
spectrum for the coupled system with a gap of magni- 
tude 2A. Using perturbation theory, they showed for 
A > Uk that Ha ^ 2/AJ2k'^kaiak = {2/A)Ha and 
-Hb ^ 2/AY.^{~iJk)blbk = {2/A)Hb. The thermal cor- 
relations obtained here can be thus regarded as direct 
consequence of this result when we apply it to a quan- 
tum quench where the coupling A is switched off at t = 
and we exploit the relationship between the GGE and the 
reduced density matrices Pa(b) described in section ID 



Another interesting consequence of the above result 
is the possibility to use quantum quenches to prepare 
systems with exactly solvable dynamics in states whose 
correlations will become indistinguishable from thermal 
correlations after the quench. However, it is unfortunate 
that the requirement of a large gap (i.e. the condition 
that A ^ 1), implies that thermal states that can be thus 
obtained are characterized by extremely high effective 
temperatures (cf. Eq. 



III. COUNTER-EXAMPLE 

The above result on the emergence of thermal behavior 
at long times stems from the existence of strong bi-partite 
entanglement in the initial state. However, as we show 
in this section, the existence of such entanglement is not 
a sufficient condition for the emergence of thermal corre- 
lations. Indeed, whenever the initial state is gapless, no 
thermal behavior can be expected even for large entan- 
glement. The Luttinger model [29[ is one such example, 
as we show below. 

Let us consider a quantum quench in the Luttinger 
model (LM) [l3|. The initial state is assume to be a 
mixed thermal state po = Z^^e'^" 1^ . The initial 



Hamiltonian H, 



LM 



H + i/int , where 



H 



LM 



2'KVp 



J2PR{k)PR{~k) + pL{-k)pL{k), (38) 



A;>0 



i^int = -J2PRWPLi-k) + PR{-k)pL{k). (39) 



L 

k>0 

In the above expression, pmL) is the density of the 
right (left) moving fermions in the Luttinger model, 
which propagate with Fermi velocity {vp) [lO[- The 
densities obey the commutation rule [pa{k), pp(k')\ = 
A:L/27r(Sfc-|_fc'(5a,^(a, /? = R,L). Therefore we can define 
two pairs of bosonic operators: PL{k) ~ \JkL/2TTa\, 
PR{k) = y^kL/2Trbl and PL{-k) = p^ik), PR{-k) = 
p\{{k). Thus, the Hamiltonian describing the system at 
t < can be written as 



H, 



LM 



E 

fc>0 



vpk^alak + bibk) + T^-kialbl + Ukbk) 

ZTT 



(40) 

Note that the LM is different from the example that has 
been discussed in sectionllllbecause the term that couples 
the two subsystems is proportional to fc, while in the 
previous example (cf. Eq. llSp it was a constant. This 
makes the Hamiltonian Hq'^^ gapless, as we show in the 
following paragraph. 

A standard way to diagonalize (|40)) is to introduce 
a bosonic canonical transformation: Ak = cosh (f>k0.k — 
sinh(/)fc&|, and Sk = — smh(f>kak + cosh(/)fc&J,. Choosing 
tanh {2(j}k) = — A/27r. The initial Hamiltonian reads now 



H^'' = Y,nk{AiAk + BlBk), 



(41) 



fc>0 



where ^k — WFfc-\/(-'- ~ {A/2ttvp)^). Thus, the energy 
spectrum of Hq^' is is gapless. 

According to the discussion of section |T] and using the 
methods of Ref. [l^l , after turning off the interaction de- 
scribed by iJjnt oc A at i = 0, the asymptotic behavior 
of the correlations can be described by a GGE matrix, 
which can be written as pgge as in Eq. ^ with a{k) and 
/3(fc) given by the entanglement spectrum of the subsys- 
tems A and _B, of modes a^, a\ and 6^, h\, respectively. 



The occupation numbers are 
sinh^(0fc) = ] 
in Eq. ^ are: 



^''(fc) 



'(fc) 



i/2{vFk/Q.k - 1). Hence, aik 


) and /3(fc) 


"(fc) = In — i FT = ^' 

yvpk- ilk J 


(42) 


R(i\ 1 (vFk + nk\ 
P{k) = In — - — — = e, 
\VFk-ilkJ 


(43) 



where e = 2[ln(27rfF/A + ^J{2t:vf/AY ~ 1)] is a con- 
stant. Hence, the entanglement Hamiltonian take the 
form: Ha = Z^fc^^^'^fc ^nd Hb = Z^fc^^^^^- Thus, 
we find that 'Ha{b) is not proportional to Ha(b)- It 
then follows that he density matrix of GGE, pgge = 
^GGE e"^^-*"*"^^^ does no longer reduce to a thermal 
ensemble. Thus, we conclude that the existence of bi- 
partite entanglement in the initial state is not a sufficient 
condition for the emergence of asymptotic thermal cor- 
relations. An additional condition, such us the existence 
of a gap, is appears to be required. 



IV. DISCUSSION AND CONCLUSION 

In this section, we would like to discuss a number of 
points concerning the above results. The first point is 
concerns the opposite situation to the one discussed in 
section [TTl We could ask what happens if initially the 
two subsystems A and B are not coupled and, at i = 0, 
they suddenly become coupled so that entanglement is 
created. Can we still expect the asymptotic correlations 
following such a quench to be described by an essentially 
thermal ensemble in a certain parameter regime? The 
answer is no, as we explain below. 

To address the above question, let us imagine that, ini- 
tially, the bosons are free to hop everywhere and there 
is no super lattice. However, at t = a superlattice is 
imposed, say by the sudden application of an extra pair 
of counter-propagating laser beams. Mathematically, the 
system at i < is described by the Hamiltonian H (cf. 
Eq. [501) ^'nd its subsequent evolution at i > is described 
by Hq (cf. Eq.IH]). The density matrix of GGE is deter- 
mined by the occupation numbers rfik) = Trpo7fe7fe = 
sin^6'fc and n^{k) = Tr poSlSk = cos^6'fe, for a half-filled 
lattice. Proceeding as above, the GGE density matrix 
describing the asymptotic correlations is: 



Pgge 



Z. 



GGE 



exp 



E 



^hh^-sls.) 



(44) 



for A ^ 1. Although the above result may appear to be 
a thermal ensemble, we must recall the dispersion of the 
eigenmodes is not a;^ = — 2cosfc but E^ = ^/uJl + A^. 
Thus, the temperature becomes again mode dependent 
and equal to T{k) = AEk/uJk, which is consistent with 
the numerical results of Ref. [y] , where lack of thermaliza- 
tion to the standard Gibbs ensemble but thermalization 
to the GGE was numerically found in this case. 

Thus, it appears again that the emergence of thermal 
behavior in exactly solvable models requires, at least in 
the simplest case, the existence of an energy gap in the 
spectrum of the Hamiltonian that determines the initial 
state. However, as we have already briefly mentioned 
in section [Hi the thermal ensemble is just an approxi- 
mation to the more general GGE, which applies in all 
circumstances. Indeed, the GGE can reproduce the be- 
havior of the asymptotic correlation functions but it can- 
not reproduce the behavior of all observables [lOl . This is 
because, in its simplest version of Eq. (O, the GGE does 
not capture all the correlations between the eigenmodes 
that exist in the initial state. For example, in the exam- 
ple of sectionini {Iklk+T,)GGE = {IklGGE X {Ik+7r)GGE 7^ 

{^olIklk+Trl^o), where |(f>o) is the initial state (i.e. the 

ground state of Ho, Eq.[Tl]) and I k+K = fl+xf^+K- This 
has important consequences, for example, when consid- 
ering the energy fluctuations: 



a^ = ($o|i^'|*o> 



- {^o\H\^o? 



.2 V^ ^fc 



(45) 
(46) 



On the other hand, if we compute the same quantity 
using the GGE, we find: 



<^GGE-^ 2^2^ 



1 



(47) 



As we have shown in section |TT1 for A ^ 1, the GGE 
tends to a thermal Gibbs ensemble (TGE) with Toff = 
A/2 = I3~^ {pGGE -^ Ptge)- And according to the 
fluctuation-dissipation theorem, for a thermal ensemble 
Ptge at a temperature Teff , 



Cv 



2 d^h\ZTGE 



dPf 



eff 



V,... 






(48) 



where Cy the heat capacity of the system. Yet, the ac- 
tual energy fluctuations of the system following a quan- 
tum quench in which a superlattice of strength A 3> 1 is 
turned off at i = are given by cr^ = 2aGGE — '^o'tge 
(cf. Ea.l47)). Therefore, we conclude that the fluctuation- 
dissipation theorem breaks down in the model of sec- 
tion ini in spite that the asymptotic correlations appear 
to be essentially thermal for A ^ 1. 

Nevertheless, although the effective temperature ob- 
tained in Eq. (j33p has no thermodynamic meaning in 
the sense of the fluctuation-dissipation theorem, it still 



represents a quantity that is worth determining experi- 
mentally. The reason is that T^g is measure of the en- 
tanglement in the system. Indeed, determination of Tofi 
from, say, a measurement of the boson momentum distri- 
bution, should allow for a determination of the entangle- 
ment spectrum of 'Ha(b)^ which, according to the discus- 
sion in sections |T] and |lTl is given by Ti. 



A(B) — Ha{B)/TcS- 

Thus, an experimental determination of the von Neu- 
mann entropy, Sa{b) = - Trp^(B) logj pa(b), {pa{b) = 
e~'^^'^ / Z A(B)) would be also possible. Similar remarks 
are applicable to the counter-example discussed in sec- 
tion lllli provided we exchange the role of Tes by e. In this 
case, however, we cannot expect thermal correlations. 

In summary, we have presented a simple instance of 
a quantum quench in which a quantum quench in an 
exactly solvable system can produce essentially thermal 
correlations. The emergence of thermal correlations from 
the generalized Gibbs ensemble (GGE) has been related 
to the existence of bi-partite eigenmode entanglement 
and a gap in the spectrum of the Hamiltonian that de- 
scribes the initial state. In this regard, we have also dis- 
cussed a counter-example demonstrating that thermal- 
ization does not happen if the initial state is described 
by a gapless Hamiltonian. Our results allow to establish a 
link between the GGE and the entanglement spectrum in 
exactly solvable systems with bi-partite entanglement of 
the eigenmodes. Thus, it makes it possible an experimen- 
tal measurement of the entanglement spectrum and other 
quantities derived from it (such like the von Neumann en- 
tropy), provided the asymptotic correlation functions of 
the system following a quantum quench can be measured 
experimentally. We have argued that this task becomes 
particularly simple when the GGE reduces to a thermal 
ensemble (or, when the entanglement spectrum has a rel- 
atively simple form, as in the counter-example discussed 
in section |TTT| . Finally, wc have also shown that, even if 
correlations may become essentially thermal, other quan- 
tities, such as the energy fluctuations, are not. This is 
akin to a breakdown of the fluctuation-dissipation theo- 
rem. 
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Appendix A: Quantum quenches in the sine-Gordon 
model 



The sine-Gordon model has been introduced in the 
main text [see Eq. ([M)) ] . This model has an exactly solv- 
able point ai K = 1 (the so-called Luther-Emery point), 
at which the Hamiltonian of Eq. ([M| can be conveniently 
represented as a quadratic form of fermion fields. To this 
end, we must use the following bosonization identity: 



i>a{x) 



1 



JSalpaix) 



(Al) 



where we have introduced the index Sq = 1 for a ^ R and 
Sq = — 1 for a = L, and the chiral bosonic fields 0q = 
K^^/^(f) + SaK^/^0, and iJjr^l are destruction operators 
for spinless fermions moving to the right and to the left, 
respectively. Using the above identity and, after a Fourier 
transformation, the Hamiltonian becomes: 

p 

+ AJ2 [^1{p)Mp) + ^iip)Mp)] , (A2) 
p 

with linear dispersion e^ = vp. 

Let us consider a quench in which the system is ini- 
tially prepared in the gapped ground state of Hsq, or 
more generally, in a thermal state defined by a density 
operator p = e~^'°/^ with temperature T. We then as- 
sume that the coupling g is suddenly turned off at i = 0, 
and therefore, for t > the time evolution is governed by 
Hq. In the long times regime, the expectation value and 
correlations of a broad class of operators for long times 
can be described by the GGE [31| . 

The sine-Gordon Hamiltonian at the Luther-Emery 
point can be diagonalized by the Bogoliubov transfor- 
mation 

iI^r[p) = cos 9pipc{p) - sin 9p%p^ {p) (A3) 

ipL (p) = sin 9p4'c {p) + cos Op^l^^ {p) ( A4) 



provided we choose 



A 
tan26'„ = — r. 



(A5) 



In terms of the new variables 'il)v,c it turns out to take the 
diagonal form 

H,G ^Y.^p H (P) V'c {p) - i^l {p)^v ip)] ■ (A6) 



with dispersion Ep = \/[£p]'^ + A^. This Hamiltonian is a 

continuum version of the superlattice model discussed in 
section ini which is obtained from (J19p in the limit where 
k =p —?' 0. 



are: 



: an initial thermal state, the cigenmodc occupati 
{ijl{p)Mp)) = I (l + cos20ptanh|^) (A8) 



:ions 



(A7) 



2TJ 



and hence the values of a and P that determine the GGE 
read: 



a{p) = log 
Pip) = log 



1 — cos 26p tanh ^ 
1 + cos 26p tanh ^ 

1 + cos 26p tanh ^ 



1- 



cos 29p tanh ^ 



(A9) 
(AlO) 
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